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We develop the auxiliary held formalism corresponding to a dilute system of spin- 1/2 fermions. 
This theory represents the Fermi counterpart of the BEC theory developed recently by F. Cooper 
et al. [Phys. Rev. Lett. 105, 240402 (2010)] to describe a dilute gas of Bose particles. Assuming 
tunable interactions, this formalism is appropriate for the study of the crossover from the regime 
of Bardeen-Cooper-Schriffer (BCS) pairing to the regime of Bose-Einstein condensation (BEC) in 
ultracold fermionic atom gases. We show that when applied to the Fermi case at zero temperature, 
the leading-order auxiliary field (LOAF) approximation gives the same equations as those obtained 
in the standard BCS variational picture. At finite temperature, LOAF leads to the theory discussed 
by by Sa de Melo, Randeria, and Engelbrecht [Phys. Rev. Lett. 71, 3202(1993); Phys. Rev. B 
55, 15153(1997)]. As such, LOAF provides a unified framework to study the interacting Fermi 
gas. The mean-field results discussed here can be systematically improved upon by calculating the 
one-particle irreducible (1-PI) action corrections, order by order. 
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I. INTRODUCTION 

One of the most remarkable achievements of the past 
decade concerns the advent of experimental studies in 
ultracold fermionic atom gases of the crossover from the 
regime of Bardeen-Cooper-Schriffer (BCS) weakly bound 
Cooper pairs to the regime of Bose-Einstein condensa- 
tion (BEC) of diatomic molecules [1-11]. In turn, these 
studies made possible for the first time the experimen- 
tal study of the ground-state properties of a many-body 
system composed of spin- 1/2 fermions interacting via a 
zero-range, infinite scattering length contact interaction. 
This regime is known as the unitarity limit [12] and is 
of particular interest in astrophysics because of its impli- 
cations regarding the equation of state for neutron mat- 
ter [13], thus emphasizes the far-reaching implications of 
these recent studies. 

Much of the theoretical work on systems composed of 
spin-1/2 fermions interacting via an adjustable, attrac- 
tive potential has focussed on interactions that are gov- 
erned by a single parameter, namely the s-wave scatter- 
ing length, a, of two atoms with different spin compo- 
nents [14, 15]. This description is valid only if \a\ ro 
and kpr <C 1, where r is the range of the potential 
and kp denotes the Fermi momentum of the gas and is 
conventionally related to the total density of particles, p, 
by the noninteracting Fermi gas formula 
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The above momentum integral is performed over the in- 
terior volume of the Fermi sphere and a denotes the 
spin component of the fermion, i.e. a = ±1/2. Then, 
the only independent dimcnsionlcss variable in the prob- 
lem is akp. Thus, this description is only applicable to 
dilute systems like the ultracold fermionic atom gases, 
not the high-density regime found in conventional su- 
perconductors [16]. The potential supports a 2-body 
bound state for (akp) -1 > 0, but this molecular state 
passes through zero energy and vanishes into the con- 
tinuum at (akp)^ 1 — 0, the position of the Fesh- 
bach resonance. Then, the BCS and BEC limits cor- 
respond to (a kp)^ 1 — >— oo and (a /cp) -1 — >-+oo, respec- 
tively, whereas the unitarity limit is defined as the limit 
near Feshbach resonances where a is much larger than the 
inter-particle distance (|a| fcp S> 1) and corresponds to 
the BCS to BEC crossover at the singularity of the scat- 
tering length. This limit is the same when approached 
with positive or negative scattering length. In the uni- 
tarity limit, the correlations are deemed to be significant 
and the system attains universal behavior, independent 
of the shape of the potential and dependent only on the 
particle density and the system dimensionality [17-19]. 

The importance of correlations in the ground state of 
dilute fermionic matter in the unitarity limit is measured 
by the numerical value of the ratio, e/sq, where e and 
So denote the ground-state energies per particle of the 
interacting and noninteracting systems, respectively. We 
recall the noninteracing energy density is defined as 




where we introduced the notation, = jk 2 with 7 = 
h 2 /(2m), and £f = 7fcp ^ s * ne Fermi energy. 

An upper bound to the value of e/eo at zero tempera- 
ture was set by the quantum Monte Carlo (QMC) study 
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performed by Carlson et al. [20-23] that gave the value 
(e/eo)QMC = 0.40(01). Instead, the "universal" curve de- 
scribing the BCS to BEC crossover in the standard BCS 
variational picture [24, 25] derived by Leggett [14] gives 
the BCS value of the chemical potential to the Fermi en- 
ergy ratio, ((J-/£f)bcs = 0.59. Other theoretical and ex- 
perimental values for e/eq are summarized elsewhere [26- 
29]. 

Recently we introduced a new theoretical framework 
for a dilute gas of Bose particles with tunable inter- 
actions [30], the bosonic counterpart of the system of 
fermions discussed in this paper. For the Bose system, 
our theoretical description is based on a loop expansion 
of the one-particle irreducible (1-PI) effective action in 
terms of composite-field propagators by rewriting the La- 
grangian in terms of auxiliary fields related to the normal 
and anomalous densities [30] . The leading-order auxiliary 
field (LOAF) approximation in the case of an interact- 
ing dilute Bose gas describes a large interval of values of 
the coupling constant, satisfies Goldstone's theorem and 
yields a Bose-Einstein transition that is second order, 
while also predicting reasonable values for the depletion. 

In this paper we will derive the corresponding auxil- 
iary field formalism for a dilute fermionic atom gas with 
tunable interactions, thus establishing the generality of 
our auxiliary field formalism. We will show that the 
LOAF approximation in the fermionic case corresponds 
to the BCS ansatz. At zero temperature the fermionic 
LOAF equations are the same as the equations derived 
by Leggett [14], whereas the finite-temperature results 
correspond to those discussed earlier by Sa de Melo, Ran- 
deria, and Engelbrecht [24, 31]. Hence, we find that the 
BCS ansatz is the only relevant auxiliary field theory in 
a dilute interacting Fermi gas. In the auxiliary field ap- 
proach, one can systematically improve upon the LOAF 
approximation by calculating the 1-PI action corrections, 
order by order. A related approach for the relativistic 
four-fermi theory can be found in Refs. 32-34. 

This paper is organized as follows: In Sec. II we dis- 
cuss the partition function for an infinite homogeneous 
system of spin- 1/2 fermions with arbitrary populations 
of spin-up and spin-down fermions. In Sec. Ill we dis- 
cuss rewriting the Lagrangian in terms of auxiliary fields. 
The resulting effective action is discussed in Sec. IV. The 
corresponding properties of an uniform system in equi- 
librium are derived in Sec. V. In Sec. VI, we specialize 
to the case of systems with equal populations of spin-up 
and spin-down fermions. We summarize our findings in 
Sec. VII. 



II. THE PARTITION FUNCTION AND PATH 
INTEGRALS 

For a grand canonical ensemble, the partition function 
for a collection of interacting Fermi particles can be writ- 
ten as 



where Cl[T,fi,V] is the grand potential and we have set 
ft = 1/T, with the temperature measured in fee units. 
Here T is the temperature, /i the chemical potential, and 
V the volume. Using the second law of thermodynamics, 
we find 



N[T,ft,V] = - 

p[T,ft,V\ = - 

S[T,n,V] = - 
and the energy E is given by 



/dQ\ 
V d/j, Jt,v ' 



(2) 



(3) 



E = n + TS + (iN . 
The partition function can be written as a path integral, 

Z[T, (jt,V]=ftf [ [ DV'D^* e S[^r;T,^v] ; (4 ) 



where S[ip,ip* ;T, [i,V] the negative of the Euclidian 
(thermal) action obtained by mapping the physical ac- 
tion to imaginary time, 1 1— > —ihr. We consider here the 
action S for a collection of fermions interacting by means 
of a short-range contact potential V(r) = Ao <5(r) is given 

by 



(5) 



where we have put fd 3 x f Q dr, and where 



(6) 



We have suppressed the dependence of quantities on 
the thermodynamic variables (T, H±,V). The fields 
are described by two-component complex anticommut- 
ing Grassmann fields ip a (x) which obey the algebra, 

Mx),rA x ')} = M x )^A x ')} (7) 
= {C(*),C'(*')} = o, 

with x = (r, tj, and a = ±1 correspond to the usual 
spin-up (f) and spin-down (!) fermions. Using a left 
derivative convention for Grassmann derivatives, varia- 
tion of the action with respect to ip^(x) leads to thermal 
equations of motion, 



{-7V 2 



Z[T,fx,V]=e 



-PQ[T,im,V] 



Tr[ e 



-0 ( H—fi N ) ' 



(1) 



— - n a + Ao [i}-„{x) ifj- a (x) ] | ip a (x) = . 

(No sum over a here.) Particle densities p a {x) and K a {x) 
are defined by 

PA X ) = CO) 1><j( x ) = P*a( x ) > ( 8a ) 

K<r(x) = 4> a {x) Ip-vix) = -K- a (x) . (8b) 
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which have the property that 
Y P<r( x ) P-A x ) = Y P*A X ) P-A x ) 



and 



(9a) 



= +Y ' 

^K a (x)K*_ cr (x)=Y2 K >Z{ x ) K -*( X ) ( 9b ) 
(7 a 

= -Y j rA x W-A x )^-A x )^A x )- 

We define k(x) = n+(x) = —K-{x). The densities p±(x) 
are real and independent whereas k(x) is complex. Intro- 
ducing four component basis vectors ip a {x) and ipa(x), 

= {^ + {x),4>_{x),r + {x),r-{x)) , (io) 

i>a{ x ) = (tp+(x),iljl(x),1p + (x),1p_(x)) , 

the density matrix can then be written as 

Pa{x) = lf, a (x)l/)\x) (11) 

/p+{x) o o -K*(xy 

p-(x) k*(x) 
k(x) —p+(x) 
\-k(x) -p-{x)j 



III. AUXILIARY FIELDS 

Following the Bose case [30], we use the Hubbard- 
Stratonovitch transformation [35] to introduce auxiliary 
fields for each of the densities described above in order 
to eliminate the quadratic interaction term in the La- 
grangian in favor of cubic interactions between the Fermi 
field and the auxiliary fields. There are six independent 
auxiliary fields possible in our case, two real fields x<r(x), 
and two complex fields A a (x), corresponding to the den- 
sities Pa(x) and Ka-(x) respectively. The auxiliary La- 
grangian is defined by 

£ aux [V>, x, A] = -C x ty>, X ] + £a[V>, A] , (12) 

where 

C x [ip,x] = Y[Xa( x ) - x oPa{ x ) sinO] (13a) 
x [x-a(x) - X p^ (7 (x) sm6] 
Y Yx°( x )x-v( x ) 



2A 



sm( 

Ao sin 2 9 



Yl[Xa( x ) P-v{ x ) X-<t( x )] 

(13b) 



£ A [^,A] = -L V[A CT (x)- \ K a (x) cose] (13c) 
2A ^ 

x [A*_ tr (x)-X K*_ <T (x) COS0] 



2A ; 

cose 

2~ 

a 

Aq cos 2 



Here we have introduced an angle 9. The auxiliary fields 
obey the same properties as the corresponding densities, 
so we define A(x) = A + (x) = — A_(x). So adding £ aux 
to C given in Eq. (6), eliminates the four-point Fermi 
interaction. Using the basis vectors given in Eq. (10), 
the action can be written in a compact way as 



[d.7 



[dx] [dx'] M x ) g- la b [<j>}(x, x') r(x') 



{ ( j )i {x} + lH )g i j {6) (cp(x)+p,i) 
2 A 



(14) 



+ M x )g a b3 b ( x ) + M x )s l (x) 

where the inverse Green function is given by 

g- la b [<f>](x,x') = 5(x,x') [Gv la M+V a b mx)] (15) 
with 



So 



bm = 



fh+ 
h+ 

-h- 
V 



where we have defined h± as the operators 

d 



h± = -~fV z ± 



dr ' 



(16) 



(17) 



and 



( X-{x) 

x' + {x) A'(x) 

A'* (:r) - X '_{x) 

\-A'*(x) 



-A'( x y 




-x'Ax), 



(18) 

Here we have redefined the anomalous fields by setting 

X±( x ) = X±(x) sine - p T , (19) 
A'(x) = A + (x) cose = -A_ (a;) cos6» , 

and introduced four component anomalous fields <fi l (x), a 
constant vector p % , and currents S l (x) by the definitions 

<P\x) = (x' + (x),x'.(x),A , (x),A , *(x)) , (20) 

Mo = ( P-,P+,0,0) , 
S\x)= (s + (x),s_(x),S(x),S*(x)) , 
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with 

= (*!-(x),x! + (x),A'*(x),&!{x)) , (21) 

(ioi = ( ^+,^_,0,0) , 

Si(x) = (s^(x),s + (x),S*(x),S(x)) . 

The tensors g a b and g 1 j{9) are defined by 

g\ =diag(l, 1,-1,-1), (22) 

Sf' 3 - (6») = diag( 1/ sin 2 (9, 1/ sin 2 0, 1/ cos 2 0, 1 / cos 2 6 ) . 

The Grassmann fields ijj a ( x ) and tl> a (x) are defined as 
in Eq. (10). It will be useful for notational purposes to 
also define ten component fields and currents using Greek 
indices as 

* a (x) = (r(x),<l> i (x)), (23) 

J a (x) = (j a (x),S i (x)). 

Note that the ip a (x) fields are Grassmann fields whereas 
the <t> l (x) fields are commuting fields, so the vectors 
^> a (x) and J a (x) are superfields. 

Setting x' — and then 9 — so that only A survives, 
recovers the first order Sa de Melo-Randeria-Engelbrecht 
theory[31]. 

IV. EFFECTIVE ACTION 

The partition function is now given by a path integral 
over all fields, 

Z[ J] = e - pn[J] (24) 
= Af //n D ^° D ^//lI D ^ D ^ e " S[ * ,J1 ' 

a i 

where S[ J] is given in Eq. (14) and again we have sup- 
pressed the dependence of Z, ft, and S on the thermody- 
namic variables. The thermodynamic partition function 
is found by setting the currents to zero. Thermal average 
values of the fields are given by 



{r{x)) = 

{<t> i {x)) = 



-1 SZ 



Z Sj a (x) 
-1 SZ 



Z SSi(x) 



= _sa 

j=S=0 Sj a (x) 

= JO_ 

3=S=0 SSi(x) 



3=S=0 



, (25) 



3=S=Q 



Since the action is now quadratic in the fields i/j a (x), we 
can integrate these out and obtain an effective action, 



-S Bft [<t>,J] 



(26) 



where the effective action S e s[4>, J] is given by 



SeS [ 0, J] = - 



[da 



dx] [dx']j a (x)g a b g b c [4>](x,x')f(x' 

d^TrfH^- 1 ^]^^)]] 
(4>i(x) + fi 0i )g l j (<^(x) 



2 An 



+ 4> l {x)S\x) 
(27) 



We evaluate the remaining path integral by expanding 
the effective action about a point 



S e s[4>, J ] = <5>cff[0O, ^ ] 

SS e $[(f>, J] 



(28) 



[dx] 



IJJ [dx] [d. 



, 5 2 S c s[(/), J 



(5<^(x) 5ft {x') 



4>o 



x (0 < (a;)-^(a:))(^(ir / )-0S(x / ))+." 

and evaluating the path integral by the method of steep- 
est descent. The vanishing of the first derivatives define 
the saddle point ft Q , which gives 



A 



= Po[<l>o, J](x) 



(29) 



1 



Here we have used 

[dx'] g- la b [cj>} (x, x') g b c l<f>] (x', x") = 5 a c <5(x - x") 



so that 

^ a e [0](x!,X 4 ) 

5ft(x) 



[dx 2 ] [dx 3 ] s° i e 1 oW(n,i 2 



x TXT^ ^ e «P (353)^4) 

d0i(x) 

= -.9% [4>] (xi , x) V 1 c d G d e [0] (x, x 4 ) , (30) 
and defined the constant matrices V la b by 

sg- la b [4>](x 2l xz) 



5ft{x) 



= V ia b S{x 2 -x 3 )S(x-x 2 ). (31) 



The densities Po[0o, J](x) are given by the equation, 

pi[<f> ,J](x) = ^ 0a tto,J}(x)V a b iP b {<p 0l J}(x), (32) 
with Grassmann fields ipg [j, 4>o] (x) given by, 

^[0 o ,J](x)= f[dx'}g\[ft J ](x,x')j b (x'). (33) 



The densities p Q [ft, J](x) and fields tpQ [ft, J](x) are 
functionals of both ft (x) and all the currents J a {x). We 
define the fluctuation inverse Green function 6x6 matrix 
T>~^[ft](x, x') by the second-order derivatives evaluated 
at the stationary points, 



V^[ft](x,x r ) 



S 2 S cS [<t>] 
Sft(x) Sft(x') 



(34) 



-f^(x,x') + ir,[0 o ](x,x'), 

^0 
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where the polarization tensor IL lJ [<po\(x, x') is given by 

FF [fa] (x, x') = \ Tr [ Vg [0o] (a;, a 7 ) V* [#>] (a/, z) ] . 

(35) 

So inserting these results into the path integral (26) and 
integrating over the (j? fields, the grand potential function 
is given by 



0SI[J] 

= f3fl + 5 eff [0o, j] - \tt[ ln[ V- 1 [<j> ] (x,x)\ 



(36) 



Ptoa-\JJ ^ ^ 3a(x) g\ G b c [M(x, x) f(x') 



[dx 



2A 



+ <l> i(x)S i (x) 



- - Tr [ ln[ g- 1 [0„] ] ] + ^Tr[ ln[ 2?- 1 [0o] (*,*)]] + ••• 

where 57o is an integration constant to be determined. 
Instead of writing the thermodynamic potential in terms 
of the currents J a (x) 1 we can write them in terms of 
fields ^> a (x) by Legendre transforming to the thermal 
vertex potential T [ ^ ] , 

(37) 

= p n[ j] - J [dx] { Mx) g\ j\x) + Mx) s l (x) } 

= (3n + ^JJ [dx] {dx'} Mx) G^Wbix, x') V> V) 
(4> i (x)+ IM )g i j (^(x)+^) 



[dx] 



2A 
1 



Tr[ In [V- 1 [cj>][ 



which is the classical action plus the trace-log terms. 
Currents are given by functional derivatives of T[ip,(j)\ 
with respect to the fields, 



ST[¥] 

Slpjx) 



= -g a b j\x), 



5<t>i(x) 



-S\x) 



(38) 



So the derivatives of T[ ^ ] with respect to the fields van- 
ish for zero currents. As shown in Ref. 36, the last term 
in Eq. (37) is of second order in a loop expansion of the 
effective action in terms of </>-propagators and will be ig- 
nored here. For the original path integral of Eq. (4), the 
loop expansion in terms of ip propagators is obtained by 
realizing that S is measured in units of h and one can 
evaluate the path integral as h — > by saddle point (or 
Laplace's method). The expansion in h leads to the loop 
expansion. Similarly one can insert an artificial small 
parameter e into the effective action by replacing S e g by 
Scs/e in Eq. (26). Powers of e then counts powers of 
loops in the composite field cf> propagators. After orga- 
nizing the series in e to some specified order one then sets 

€ = 1. 



V. UNIFORM SYSTEM IN THERMAL 
EQUILIBRIUM 

For the case of a uniform sample and thermal equilib- 
rium, the fields ip a and <j) 1 are independent of x = (r, r). 
In addition since the Green functions are periodic or anti- 
periodic in r, we can expand them in a Fourier series, 



GMx,x / ) = j j J2Sm^n)e^ 



(r— r )— w»(t— T )] 



(39) 



where u> n — (2n + l)//3 are the Fermi Matsubara frequen- 
cies. So using /[dx] = fiV, at thermal equilibrium and 
for uniform systems, the thermal effective potential from 
Eq. (37) is given by 



v eff [*] = r[ *]/y = v + ^ a v[<t>] a b i> b 

2A 

Here the matrix (? _1 [<£](k, n) is given by 

G-^lflc n)-f A(k ' n) B ( k > n ) 



(40) 



with 



A(k,n) 
-A* (It, n) 

B(k, n) 
-S*(k, n) 



e k + x'- - iu n 

£fc + X+ ~ iuJ r 

~ X- ~ iuJ n 

-e k -x'+- 



-A' 
A' 



A'* 
A'* 



(41) 
(42a) 

)• 

(42b) 
(42c) 

(42d) 



Note that the matrices A and B satisfy the integration 
conditions for the Grassmann integral. From (18), the 
matrix V[<j>] is given by 



(43) 



V[0] - 


( X '~ 







x'+ A' 
A'* -x'- 


-A 






V-A' 


* 


-x' 


Then we find 










X-P- 


f + X+P- + 


A'*K 


Also we find 














2A 










-){x'-+P-] 


- + A ' 




A 


sin 2 9 





(44) 



(45) 



A 



l\2 



Ao cos 2 9 
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The grand potential per unit volume is the value of the 
effective potential evaluated at zero current, or when 



with 



= 0, 



(46) 



for all values of a. From Eq. (44), we find 

\ i> a V\<$>\\ ^ = X '-P+ + x'+P- + A' k* + A' *k , (47) 
and from Eq. (45) , we get 



2A 

(x'++M+)(x'-+/0 , |Af 



(48) 



Aq sin 



Aq cos 2 6 



In order to compute Tr[ln[£/ _1 [(/)](k, n) ] }, it is simpler 
to interchange rows and columns of the CJ -1 [0] matrix so 
as to bring it into block diagonal form. To do this, we 
redefine the fields -ip a (x) in the following way: 

r (x) h. ( i, + (x) ,ip_(x), r+ (x) , V- (x) ) , (49a) 
ifv, (as) h> ( V* (a?) ,1>-(x),il> + (x) , C (a;) ) • (49b) 

Then from Eqs. (41) and (42), the Fourier transform of 
the inverse C/-matrix is of the form, 



where 
G-\k,n) = 



G _1 (k,n) 

-G-^(k,n)J ' 



(50) 



(51) 



(52) 



e k + x'- - iu n -A' 

-A'* -e k -x'+-iu 7 

After some algebra, we find 

det[G _1 (k, n) ]=—(«! + ia 2 ) , 
det[G _1 ^(k,n)] = — (a± — ia<i ) , 

where 

«i = (efc + X+){tk + X-) + |A'| 2 + u? n , (53a) 
aa = (x / --X+)wn. (53b) 

So we find 

»[ln[0- 1 (k,n)]] =ln[det[g- 1 (k,n)]] = (54) 

= ln[det[G -1 (k,n)] det[G -1 f (k, n) ] ] 

= ln[af + al] = ln[w* + 2 6w 2 + c] 
= ln[(^ +w 2 )(cj 2 +tj 2 )] 

= ln[^+^]+ln[^+w 2 _], 
where w 2 . are the two roots of the equation 



u4 - 2 6w| + c = 



l=b±y / b 2 -c, (55) 



& = ( e , + x V)( efc + X '_) + |A'| 2 + - X-)72 (56) 



[(e fc +xV) 2 + (^+X-) 2 + 2|A'| 2 

2 



2 

c= [(6 fe + xV)(e fc +x'-) + |A'| 2 



The square root of the discriminant is given by 
\/b 2 -c (57) 
= \W+- X- 1 yfi (ek + X'+) + (e k + X '-) ] 2 + 4 | A'| 2 . 
so the frequencies w±, which depend on k, are given by 

<4 = ^ { [ + x' + ) 2 + (<* + X'J 2 + 2 | A'| 2 ] (58) 
± lx'+ - X'-I 0(e fe +x' + ) + (e fe + x'-)] 2 +4|A'| 2 } . 
So from (54) 



■TrtlnlCT^kjTi)]; 



(59) 



/c\ 3 k . 
— ^ln[det[^(k,n)]] 

/pf{EH^+^] + EH«S+^-]} 



1 /" d 3 fc 
1 /' <1 : 7> 



d 3 fc r cj + + cj_ 
(27)3 \ 2 

+ - { ln[ 1 + e"^+ ] + ln[ 1 + e~^ u - } } } . 

Then from (40), the effective potential becomes 

V cff [*] - V + x'+P- + X-P+ + A' k* + A'*k (60) 
(x'++M+)(x'- + M-) |A'| 2 



A sin 2 9 
d 3 k r cj + + cj_ 



An cos 2 



(27T) 



{ 



+ i { ln[ 1 + e-^+ ] + ln[ 1 + e^"- ] } } . 
Expanding w± in a Laurent series about k — > oo, we find 
- . 1 r-7 , .7 ! , l A T , .. (61) 



eu + \[x' + + x'-r-r — 



So using dimensional regularization [37] , the effective po- 
tential becomes 

V c ff[*]=X>-+x'-P+ + A' K *+A'* K (62) 
(XV+M+)(X-+M-) |A'| 2 



A sin 2 i 



A cos 2 



7t) 3 I 



l A '| 2 

2 efc -2 [x '+ +X '- ] ~^T 



d 3 fc c w + + cj_ 1 
(2 



{ln[l 



] + hi[ 1 + e~ Pu ~ 



!} 
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where the coupling constant is related to the s-wave scat- 
tering length, a, i.e. A = 87r7 a/. Recall that p± = ip± ip ± 
and k — ip+ ip- . 

We recover the thermodynamic grand potential by 
evaluating V e fr[ ^ ] at the minimum of the potential when 



d 1p a 



, and 



dd>i 



, (63) 



for all values of a and i. For the Grassmann ip a fields, 
derivatives of the effective potential (62) with respect to 
V»± give 

[x' + X- + |Af ]V> ± =0. (64) 
The above can be satisfied only if ip± = 0. 

VI. EQUAL CHEMICAL POTENTIALS 

In this section we set /i+ = /i_ = fi and x+ = X- = x'j 
so that only the total particle density is fixed. For this 
case the only possible solution for the Grassmann fields 
is the first case above where ip± = 0. The frequency 
spectrum is given by 



u? k =ufi=uyL = (e fe + x') 2 + |A / 



(65) 



and from (62) , the effective potential for this case is given 

by 



V„ff[tf 1 = - 



- 2 



(x' + m) 2 

A sin 2 (9 
f d 3 k r 1 

J W) 3 I 2 
1 



IA 



/|2 



/3 



ln[l 



A cos 2 9 
Wfc - £fc - x' 



(66) 



|Af 

2e fc 



]}■ 



The gap equation for the x' field is now, 



X 



_ 1 
Asm 2 (9 ~~ 2 
1 f d 3 fc 
" 2 J (2^ 



_d 3 fc 
(2' 

£fe+x' 



-1]+1 



{f*±2L[2n(^ k )-l]+l} ) 



} 

(67) 



where the non-interacting Fermi particle number factor 
n(x) is defined by 

n(a;) = l/[e x + l]. (68) 
The gap equation for A' is 

A' f d 3 k ( ( doj k \ rn , n , -i A' i 



(2 



Again, the factor of A' cancels, and we get for the A' 
gap equation 

d 3 fc 
(27p 



1 

A cos 2 (9 



3 k r 1 l-i 
_{_[2n(^ fe )-l]+-[. (70) 



The total particle density is given by 



p 



dfi 

d 3 k 



A sin 2 9 



(71) 



J_(!^L [2n(/3wfe) _ l]+1 } 

7T) 3 L UJ k i 



It is convenient to scale momenta and energies in terms of 
the Fermi momentum, fcp, and Fermi energy, e F = jk F , 
respectively. We introduce 



k = k/kF, /x = /z/e F , A = A/e F 



(72) 



Q k =u k /e F = v /(fc2+x') 2 + A2, 
T = T/e F = T/T F , /3 = e F /3 = T F /T = l/f . 
Then, the rescaled equations are 



1 



k 2 dk- 



1 l-2n{PQ- k ) 
k 2 wr 



£ COS 2 9 7T J a 

1 = £ / p d fc{l-^^Ml-2n 



X' = — £ sin 2 - // 



}, (73a) 
OS**)]}, (73b) 
(73c) 



where now w 2 , = (fc 2 + x') 2 + |A| 2 and we introduced 
£ = fcpa. Eqs. (73) are to be solved selfconsistently for \x 
and lA'l. 



A. Zero temperature (T = 0) 

At zero temperature so that n{pQ k ) — 0, Eqs. (73) 
reduce to 



1 



£ cos 2 9 ir J 



3tt 



- ' d*|l- — > , 

L Ml J 



(74) 



fc 2 dfc<^ 1 



£ sin 9 — n . 



In Fig. 1 we illustrate the solutions of the gap equations 
(74) for A' and n in reduced units as a function of l/£ 
for values of 9 = 0, 7r/12, ir/6 and 7r/4. For 6 — 0, our 
results reduce to the variational equations discussed by 
Leggett [14]. In the unitarity limit (i.e. for l/£ = 0) one 
obtains p/ep = 0.59 and A/ep — 0.69. 

For 9 0, the chemical potential has a singularity at 
l/£ = 0. This indicates that the only physical theory 
corresponds to choosing 9 = 0. Hence, the BCS theory 
is the only relevant auxiliary field theory for a dilute gas 
of fermions. 



B. Critical temperature (A' = 0) 

At finite temperature, the critical temperature and 
critical chemical potential correspond to the point where 
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FIG. 1. (Color online) Zero temperature solutions for A' and 
(j, of the gap equations in scaled units vs. l/£ for several values 
of the parameter 6 = 0. 



FIG. 2. (Color online) Solutions for T and p of the gap equa- 
tions in scaled units at the critical point (A' = 0) vs. 1 /£ for 
several values of the parameter 6 = 0. 



the gap A = 0. At the critical point, the spectrum be- 
comes Uk = £fc + x' ■ F° r this case, Eqs. (73) become 



1 



£ cos 2 9 



tanh(/3w fc /2) 



} (75) 



fc 2 dfc{l-sgn[fc 2 +x] tanh(/3wfe/2) } 



X 1 = — £ sin 2 9- p. 

Solutions of Eqs. (75) for various values of 9 are shown 
in Fig. 2. For 9 = 0, our results are the same as those 
discussed extensively by Sa de Melo, Randeria and En- 
gelbrecht in Refs. 24 and 31. For 9 ^ 0, the chemical 
potential at the critical temperature has a singularity in 
the unitarity limit, again indicating that the case of 9 = 
corresponds to the only physical theory for a dilute gas 
of fermions in the auxiliary field formalism. 



C. Thermodynamics 



From Eq. (66), the particle number density is 



N 
V 



= 2 



1 /dfl\ 
V\ dp Jt,v 

d3fc (US 



Weft 
dp 



A sin 2 9 



(76) 



where 



p(k) 



£fc + X 



[l-2n(fiuj k )\ } . 



(77) 



The zero-temperature momentum distribution of the par- 
ticle distribution function, p(k), is shown in Fig. 3 for 
several values of the parameter £. For completeness, we 
also depict the momentum dependence of the dispersion 
relations, ui^, for 9 = and the same values of the pa- 
rameter £. We note that the location of the minimum in 
the dispersion relation shifts smoothly to zero momen- 
tum and disappears for £ > £ c ssO.55, indicative of the 
crossover character of the BCS to BEC transition [25]. 

The pressure is also obtained from Eq. (66), as 



, d 3 k 

2 



-V, 



eff 



f A 6 k r 1 



(x' + fj-)' 

Asin 2 6> 

Wfc — £fc — > 

+ ~ln[l + 



v'|2 



A cos 2 9 
lA'l 2 



(78) 



2e fc 
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In scaled variables, the pressure is given by 



P 
pe F 



3tt |A 



112 



+ |jT °A 2 dfe{w s -fc 2 



(79) 



|A? 

2 k 2 



■= ln[ 1 + e 



From Eq. (2), the entropy per unit volume, s, is given by 



V 

2(3 



I -I = - = /3 2 



1 /d£l\ 

d 3 /c 
(2 
d 3 fc 



(3 2 on 
V 80 



dV eS 

3(3 



iU ( i 



/3 

{n(/3w fc ) ln[n(j8a;jfe)] 



(2tt) 3 

+ [l-n(/3w fc )] ln[l-n(/3u; fc )]} 



(80) 



or, in scaled units. 



k 2 dfc 



ws + i ln[l + e 



(81) 



From Eq. (3), the energy per unit volume, e, is given by 

e = E/V = Veff + Ts + MP (82) 
lA'l 2 



A sin 2 
d 3 fc 



A cos 2 (9 



/d. /c r 
| [ 2n(/3w fc ) - 1 ] uj k + e k + X ' 



or, in scaled units, 



e 



P 
pe F 

2 



T 



pe F e F 



(83) 



, ■ 2n 37r l A T 

— £ sm 2 9 + p - — J— ^ 
3n 8 £ cos 2 9 



>dfc{o) £ [l-2n(«] -P-x'-^f} 



Here x' and A' solutions of Eqs. (73). Comparing 
Eqs. (79) and (83), we see that at T = 0, 



-p + np. 



(84) 



For illustrative purposes in Fig. 4 we depict the zero- 
temperature pressure and energy per unit volume as a 
function of l/£, for several values of 9. The pressure 
and density have singularities in the unitarity limit, con- 
sistent with our previous results that the case of 9 = 
corresponds to the only physical theory for a dilute gas 
of fermions in the auxiliary field formalism. In Fig. 5, we 
illustrate the equation of state, E/pV, vs. l/£ for 9 = 0. 




k/k c 



FIG. 3. (Color online) Momentum dependence of the zero 
temperature particle distributions, p(k), and dispersion rela- 
tions, Uk, for 9 = and several values of the parameter l/£. 
We note that the location of the minimum in the dispersion 
relation shifts smoothly to zero momentum and disappears 
for £ > £ c ~0.55, indicative of the crossover character of the 
BCS to BEC transition. 



D. Contact interaction relations 

For fermions interacting via short-range potential, Tan 
derived a set of universal relations in Refs. 38 that are in- 
dependent of the details of the short-range interactions, 
some of which have been verified in experiments [39, 40]. 
In particular, Tan relates the fermion momentum distri- 
bution, p(k), at asymptotically large momentum to ther- 
modynamics quantities such as the energy of the system 
per unit volume: Tan showed [38] that the fermion mo- 
mentum distribution satisfies the property that 



p(k) 



C_ 
fc 4 



(85) 



in the large momentum limit, where C is the contact den- 
sity. This results was observed experimentally by Stew- 
art et al. [39]. Next, according to Tan's "adiabatic sweep" 
theorem [38], the variation of the energy per unit volume, 
e, with respect to the inverse scattering length is given 

by 



da- 1 2vr 



(86) 



This result was also verified experimentally by Stewart 
et al. [40]. 
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FIG. 4. (Color online) Zero temperature pressure and energy 
per unit volume vs. l/£ for several values of the parameter 
= 0. 



We will show here that the LOAF approximation sat- 
isfies these two Tan relations: First, from Eq. (77), we 
find that indeed 



Cloaf 

_ k^ 



p(k) 

with the LOAF contact density 
Cloaf 



(87) 



A' 2 

4^2 



Second, we take the derivative of the energy per unit 
volume, e, given in Eq. (82) with respect to the inverse 
scattering length. Using Eq. (82) and recalling that at 
the minimum we have dV e s/dx' = 0, dV c s/dA' = 0, and 
dVeff/dji — —p, we find that 



de 
da- 1 



dp 
da' 1 
A' 2 



(89) 



C 



LOAF 



_ _ 1_ 

87T7 27T 

as indicated by Tan's relation (86). 



E. Unitarity limit 

From Eqs. (73), we see that in the unitarity limit, 
l/£ — >• 0, the gap equations can only be satisfied if 6 = 0, 



> 







\ / 


\ : 







-2.0 



-1.5 



-1.0 



-0.5 0.0 



0.5 



1.0 



FIG. 5. (Color online) Zero temperature equation of state for 
6 — vs. l/£. We note that the ratio e/p is equal to 3/2 both 
in the limit of a noninteracting Fermi gas and in the unitarity 
limit. 



which shows yet again that the case of 9 — corresponds 
to the only physical theory for a dilute gas of fermions 
in the auxiliary field formalism. Hence, in the unitarity 
limit, the scaled gap equations become 

o = ifd*{i-^ 

7T Jo LUk 

i = 6 / k 2 dkh- — 

2 .In l - Wfc 



fl [l-2n(j3u k )]\, (90a) 

I u) k J 

/ v r k 2 - ii -i 

/ k 2 dk\l "\l-2n(Bu k )]\, (90b) 

Jo L w fc ' > 



where now u 2 , = (fc 2 — p) 2 +\ A| 2 , and we have dropped the 
bar notation. The pressure and energy per unit colume 
are now given by 



P 
pe F 



J k 2 dk{co k - k 2 



2 k 2 



e 3 
pe F 2 



1 ( 

k 2 dki u} k [l-2n(pu) k )} 

} 



(91a) 



(91b) 



2fc 2 



By parts integration, we have 

-3u k 1 



f 1 
/ k 2 dk - ln[ 1 

Jo P 



(92) 



2 [°° l2A1 k 2 (k 2 ~p) 
= - k 2 dk — ^ ^ n(j3u k ) . 

Substituting this into Eq. (91a), the pressure can be writ- 
ten as 



P 
pe F 



k 2 dk 



f 3 T ,2 \^ 
{- 2 [^-k +p-- 

ifiUh) } 



(93) 



2 fc2 ( fc2 ~^ 
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For the energy expression, multiply Eq. (90b) by p and 
substitute the result into Eq. (91b). This gives for the 
energy 



e _ 3 
pe F 2 



fc 2 dfc(- ^ + M(fc2 ~ M) [l-2n(/3^)- 



2fc< 



} 



(94) 



Now form the quantity 
2e - 3p 



pe F 
3 
" 2 



3 / Fdfc 



Wfe - k + /J, ■ 



JAP 
2fc 2 



2|A| 2 



Now integrate by parts, 



k dk 



Uk - k + fi - 



2 k 2 



2 k 2 
(95) 

(96) 



k dk 



., k 2 (k 2 ~p) |A| 2 



2 k 2 



Inserting this into (95) gives 



2e - ip 
pe F 



3|A|2 r dfc j 1 _* a [i-^)]\ =0 , 

(97) 

where we have used the gap equation (90a) . So this shows 
that at the unitary limit. 



2 F 



(98) 



for all temperatures T (see e.g. Ref. 41). In Fig. 5 we 
show numerically that this relation holds for T = 0. 

Using Eqs. (84) and (98), we find the unitarity limit 
results at zero temperature, 



pe F 



(99) 



in reduced units. But, at T = 0, we have p/e F = 0.59. 
Therefore, introducing the energy per particle 



E 
N 



(100) 



we obtain that at zero temperature, in the unitarity limit, 
we have 



O/^o) 



LOAF 



0.59 . 



(101) 



VII. CONCLUSIONS 



To summarize, in this paper we derived the aux- 
iliary field formalism for a dilute fermionic atom gas 
with tunable interactions. This formalism represents 
the fermionic counterpart of a similar auxiliary field for- 
malism introduced recently to describe the properties 
of a dilute gas of Bose particles [30]. Here we demon- 
strate that at zero temperature, the fermionic LOAF 
equations are the same as the equations derived by 
Leggett [14], whereas the finite-temperature results cor- 
respond to those discussed earlier by Sa de Melo, Ran- 
deria, and Engelbrecht [24, 31]. The LOAF formalism 
shows that the BCS ansatz represents the only physical 
auxiliary field theory for a dilute Fermi gas. Furthermore, 
we showed that LOAF satisfies Tan's relation regarding 
the momentum distribution of fermions at asymptoti- 
cally large momenta and Tan's "adiabatic sweep" the- 
orem. Just like in the Bose case, the auxiliary field ap- 
proach for fermions provides a systematic framework that 
allows one to improve the LOAF results by calculating 
the 1-PI action corrections, order by order. 
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